Phase Transition, Entanglement and Squeezing in a Triple- Well Condensate 
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We provide an in-depth characterization of a three modes Bose-Einstein condensate trapped 
in a symmetric circular triple well potential. We analyze how a subsystem independent measure 
of entanglement, the purity related to the su(3) algebra, scales for increasing number of atoms 
and signals correctly the quantum phase transition between two dynamical regimes in a specific 
arrangement. Moreover, this measure, which is intrinsically related to particle entanglement, also 
depicts if some squeezing is occurring when we consider the system's ground state. 
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Entanglement has played an important role for the un- 
derstanding of quantum many body aspects pQ that tra- 
ditionally belonged to statistical mechanics and quantum 
field theory. Several investigations in quantum critical 
models at T = have shown that complex entangled 
ground state contains all the important correlations that 
give rise to different phases known to exist in several sys- 
tems [2]. Thus, it is a fact that entanglement study in 
many body systems allows a deeper characterization of 
the ground state of the system undergoing a quantum 
phase transition (QPT), particularly its order. Charac- 
terization of a QPT via pairwise and collective multi- 
partite entanglement has been given in a very conclusive 
way in Refs. [3j [4], but are dependent on the specific 
partition employed, i.e., they are subsystem dependent. 
Another way of investigating subsystem independent en- 
tanglement in many-particle systems is through the gen- 
eralized purity associated to the pertinent algebra [5]. 
Beyond being a measure for the quality of semiclassical 
approximation [6] , this measure is also related to squeez- 
ing of moments of the generators of the pertinent alge- 
bra. Recently both spin squeezing and entanglement 
have been demonstrated for a 87 Rb condensate trapped 
in double and multiple wells of an optical lattice [7 , and 
subsystem independent entanglement theoretically inves- 
tigated in Ref. [6 for a double- well trapped condensate. 
In fact although those results were developed indepen- 
dently they are profoundly complementary since they 
relate QPT, entanglement and squeezing, for a system 
which is a particular realization of the Lipkin-Meshkov- 
Glick model [8, 9 . The interplay between entanglement 
and squeezing has been investigated previously in many 
instances [10 , and now seems to play an important role 
in QPT involving many bosons as well. 

In this Letter we investigate in detail a BEC of at- 
tractively interacting neutral atoms trapped in a sym- 
metric triple well potential in a three mode approxima- 
tion and show that the ground state of the model under- 
goes a QPT. A time dependent variational principle using 
SU (3) coherent state allows for a system of semiclassical 
equations that enables one to find the fixed points of 
the model and to investigate how the lowest energy fixed 
points change as the collision parameters of the model are 



varied. Since the lowest energy state in this system corre- 
sponds to a twin condensate fixed point, where effectively 
the system behaves as if composed of two wells (although 
coherences between the twin modes are still present) we 
can show that the SU (3) reduces to SU (2) while in this 
regime. The non-linear components of the Hamiltonian 
provided by the interactions of our model lead to spin 
squeezing for the ground state as we vary the scattering 
parameter, which on its turn is associated to variance 
squeezing of some of the generators of the algebra, and 
shows non-trivial phenomena such as particle entangle- 
ment related to QPT and to distinct dynamical regimes. 
This special arrangement is promising for experimental 
investigation with existent trapping technology. 

The model of a Bose-Einstein condensate trapped in a 
triple well potential has already been considered in lit- 
erature in many distinct configurations [TTJ [12]. Con- 
sidering that the coupling between the potential wells is 
sufficiently weak, we shall use the usual local modes ap- 
proximation [13] associated to the states \uj), j = 1, 2, 3, 
representing the ground states of harmonic approxima- 
tions around each minimum of the trap. Considering 
this approximation to be valid even when the total num- 
ber of trapped bosons TV is large, the overlap e = (ui\uj), 
(i 7^ j) must be very small. Keeping terms up to O(ez), 
and using the fact that the total number of trapped par- 
ticles N is conserved, we can write the Hamiltonian: 

H = ft ^2 a \ a 3 + K a T°H ~ 2 ^ ^2 <4 a i a ] a k'i (1) 

where ft = ft + 2 A (TV — 1), and the operator aj(at), 
j = 1,2,3, annihilates (creates) a boson in the state 
\uj). In the last equation we also defined the param- 
eters of the system: ft is the tunneling rate, ft is the 
effective tunneling rate, n is the self-collision parameter 
and A = <c ft is the cross-collision rate, which is pro- 
portional to the interaction frequency between bosons in 
different sites. Although A is a lower order parameter 
of the model, we see that its presence in the effective 
tunneling rate ft' = ft + 2A(7V- 1) is relevant for N > 1. 

Since the Hamiltonian (fiT) preserves A 7 ", we can take 
advantage of the homomorphism between the commu- 
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tation relations of SU(3) group generators and bilinear 
combinations of creation and annihilation operators in 
the three local modes, making use of an extension of 
Schwinger's pseudo-spin method [T4] : 



Qi = -(^1^1 — <4 a 2)> Q2 = o ( a i a i + a 2 a 2 _ 2*23(23); 



J k =i{a\dj -a\aj), 



a k a 3 



(2) 



for k = 1,2,3 and j = (k + l)mod3 + 1. Using these 
eight generators we can rewrite the Hamiltonian: 



H 



;(4Q? + 3gl) 



+A [2Qi(Pi - P 3 ) + Q 2 (2P 2 - Pi - P 3 )] , (3) 



showing that SU (3) is the pertinent dynamical group. 
Note that the tunneling term is linear in the generators, 
while the collision terms are quadratic. 

A semiclassical dynamics of the condensate can be de- 
veloped by using the time dependent variational principle 
[T5] with coherent states of the group SU (3) [16] as test 
functions, allowing us to express the semiclassical Hamil- 
tonian as 7~L(w*,w) = (TV; w\H\N; w), with the coherent 
states expanded on the Fock space basis for the three 
local modes as 



\N;w) = c 



ni+n2+n 3 =N 
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where C = 1 /(|^i | 2 + |w 2 | 2 + l)" 5 " • The vector w = 
(1^1,1^2) £ C 2 parametrizes the non-linear subspace 



of Hilbert space consisting only of coherent states and 
also represents a point in phase space associated with 
the classical analogue system [17 . The semiclassical 
equations of motion are then brought to the Hamil- 
ton's canonical form by the transformation of variables 
Wj = ^I 3 /{N-h-I 2 )e- 



1,2. While L is the 



average occupancy in the j-th local mode, <j)j represents 
the difference in phase between the macroscopic conden- 
sates located in the j-th and the third wells. The semi- 
classical approximation becomes exact if the Hamiltonian 
is linear in the generators of the group or in the classical 
limit of the system, which coincides with the macroscopic 
limit N —> 00 [18 . The phase space position and the sta- 
bility of the equilibrium points of the equations of motion 
depend only on the model parameters x = and 



A(AT-l) 



(See Ref. [19]). 



As previously noticed in Ref. [12] there are fixed point 
solutions corresponding to configurations where two lo- 
calized condensates are in phase and have the same occu- 
pational average, which are described by the three equiv- 
alent conditions W\ = u> 2 , W\ = 1 or w 2 = 1, known as 
restrictions of twin condensates. Many of the fixed points 
of the semiclassical dynamics are contained in this sub- 
regime under the additional condition of 0i,02 = 0,7r, 
which implies that w G R 2 . Without loss of generality, we 
choose only the condition w = W\ = W2 G R, so that the 
SU (3) coherent states reduce to coherent states of SU (2), 
which thus brings many features observed for a two mode 
condensate, such as the Rabi Oscillation (RO) of popula- 
tion and the macroscopic self-trapping (MST) of popula- 
tion [20] . The phase space associated with the integrable 
sub-regime of twin condensates is isomorphic to S 2 , a 
space that parametrizes the set of SU (2) coherent states 
I J = y,r = V2w = tan |e _z ^). An interesting quantity 
for the dynamical analysis in the twin condensates sub- 
regime is the population balance I z 
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Figure 1: (color online) Semiclassical trajectories in the sub-regime of twin condensates on the unit sphere for various initial conditions. 
For the three spheres we use TV = 30, Q = —1 and /i = 0, while the values of the reduced self-collision parameter are: (a) % = 1.5, (b) 
X = 1-98 and (c) % — 3. The yellow circles (red crosses) indicate points of stable (unstable) equilibrium for the full dynamics. The critical 
parameter of bifurcation in the absence of cross-collisions is X+(£t = 0) « 1.97. 
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The point w\ — W2 — 1 is a solution of the fixed point equations, independent of the values of \ an d /i, which 
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we call 1 + . Considering always /i <C x, we have an- 
other solution if x < X+(aO> which we call 2 + ; whereas 
if X > X+(aO two further solutions of the fixed point 
equations exists. In other words, a saddle-node bifurca- 
tion with critical parameter X+(aO is responsible for the 
appearance of the fixed points 3 + and 4 + . 

Fig. [I] a shows the semiclassical dynamics of twin con- 
densates in the absence of the fixed points 3 + and 4 + . We 
see that all trajectories are around the fixed point 1 + and 
therefore vary around the value I z = - , which represents 
same average occupancy in the three local modes. This is 
exactly the RO dynamical regime, where the system does 
not show any preferential mode occupation. In figure[l]b, 
where x exceeds the critical value of bifurcation, we ob- 
serve the emergence of new types of orbits in the system, 
especially around the new stable equilibrium point 4 + . 
These new trajectories oscillate around negative values of 
I z and show a suppression of tunneling between the twin 
condensates and the solitary one, resulting in the prefer- 
ential occupation of the third mode. Such a dynamical 
effect is known as macroscopic self-trapping (MST), sim- 
ilar to the case of two modes p~3j |21] . As the value of 
X increases, the region of phase space occupied by orbits 
associated with the MST regime also grows larger, as we 
see in figure [I] c. 

A quantitative error in the semiclassical approximation 
is originated when we force the state to evolve preserving 
minimal uncertainty on the space phase, i.e., obligating 
the system to remain as a coherent state. If the quan- 
tum dynamics of the system leads to an increase of the 
state uncertainty, we cannot expect a high quantitative 
accuracy in the semiclassical treatment. Including the 
situation where strong squeezing or spreading over the 
phase space occurs, leading thus to a large deviation of 
the coherent state. An effective measure of the quality of 
the semiclassical approximation is the generalized purity 
associated with the su(3) algebra [5], defined by 



P*u(3)(M) = 



9 ( WQiW 2 

N 2 \ 3 
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(5) 



This measure is derived from the total uncertainty of 
the algebra [22] and has maximum value V su (3)(\il))) = 1 
only if is a coherent state of SU(3), corresponding 
to the state of minimum total uncertainty. In contrast, 
the purity presents a minimum value ^ sn (3) (IV')) = for 
states of largest uncertainty. Therefore, V su (z) can be 
used to measure the "distance" of a particular state to 
the subspace formed only by the coherent states. That is, 
the larger the uncertainty of a state during its evolution, 
the lower is V su ^s) and less accurate is the semiclassical 
approximation using only coherent states [6]. 

In addition the generalized purity is also a genuine 
measure of separability in systems of many identical par- 
ticles. Under a basis transformation of the single-particle 
Hilbert space, every coherent state of the type Q can be 
rewritten as a separable state in each boson: 



\N;w) 



*N\ 



w\a\ 



w 2 



t \N 



1)4 

N 



N 



■|0> = |JV*> 



10} 



(6) 



where |^) 



W iK)+ W2 |h 2 )+|h 3 ) ig an arbit gtate of 

(kl| 2 + k 2 | 2 + l)2 

the single-particle Hilbert space [25] with associated an- 
nihilation operator a^. Due to the symmetrization prin- 
ciple, we see that ([6| is the only possibility of a com- 
pletely separable state in each particle, considering sys- 
tems of many identical bosons. Therefore, the states 
with V su (3)(\ip)) = 1 represent the separable states of N 
bosons, while the decrease of purity from this maximum 
value indicates the entanglement among the particles of 
the system. In this sense when the state of the system 
is less localized in the phase space, V su (^){\iji)) < 1, it 
is more entangled. Several dynamical processes can be 
responsible for the loss of purity of an initially coherent 
state. One of such processes is the squeezing of the state, 
and it is related to regular dynamical regimes of the semi- 
classical model [20]. In the semiclassical approximation 
we also observe chaos [19 in the three mode model, ac- 
companied by other processes of purity loss. 

The generalized purity can also be written as the trace 
of a squared reduced density operator, however this re- 
duction is done with respect to the algebra. Consid- 
ering p = we have V 3U (3)(\^)) = ^ Tr (pL(3))> 
where K is a normalization factor dependent on N ', 
Psu(3) = is the density operator reduced 
3 

in su(3) and {Aj} is a basis of su{3) orthonormalized in 
relation to the trace in the Fock space of N particles. 
Therefore, the entanglement can be described as the de- 
coherence of the system in its dynamical algebra. Sum- 
marizing, pure (mixed) states in the algebra are equiva- 
lent to separable (entangled) states in each boson. The 
semiclassical approximation is accurate only if the state 
remains separable in each boson, when we can completely 
describe the system by only one trajectory in the phase 
space that extremizes an action functional. Then, the en- 
tanglement is the fundamental quantum feature respon- 
sible for breaking-off the classicality of the system, and 
it signals QPT as we discuss now. 

QPT results from non-analyticities in the ground state 
energy as a function of a real parameter of the Hamilto- 
nian, characterized only in the macroscopic limit N — >• 00 
and at zero temperature. However, as stressed earlier, 
the macroscopic limit of the model is equivalent to its 
classical limit, so that the minimum of the semiclassical 
energy per particle ^, which is independent of A", can 
be used exactly to study the QPT. The extremes of T~L 
are the equilibrium points of the model. In Fig. [2] we 
show H calculated in each of its extremes as a function 
of x? with N = 30, Q = — 1 and // = 0. The occur- 
rence of a first order QPT is evident at the critical value 
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Figure 2: (color online) W, calculated at its extremes as a function 
of the self-collision parameter k showing the fixed points for the 
twin-condensate regime for TV = 30, Q = — 1 and /i = 0. 



Xc = 2, since there is a discontinuity of the first deriva- 
tive of the ground state energy, resulting from the energy 
level crossing of the fixed points 1 + and 4 + . 




Figure 3: Occupational Husimi representation Qi(I\,l2) (above) 
and phase distribution <E>(</>i , fa) (below) of the ground state for 
TV = 30, O = — 1, fi = and self-collision parameter value: (a) 
X = 0, (b) x — 2 and (c) x = 3. Remark the squeezing in phase 
while the distribution in number broadens. The fragmentation of 
the phase in (c) characterizes the QPT. 

As already shown, the fixed point 1 + is the center of 
the orbits related to the RO dynamical regime, charac- 



terized by no preferential occupation of the local modes. 
The three equivalent equilibrium points 4 + , one in each 
sub-regime of twin condensates, are responsible for the 
emergence of semiclassical trajectories associated with 
the MST regime in the solitary mode. Therefore, we 
expect the ground state to display characteristics of the 
RO (MST) regime for x < Xc (x > Xc)- Furthermore, 
the behavior of the low energy equilibrium points sug- 
gests a non-degenerate (triply degenerate) ground level 
when the self-collision parameter is lower (higher) than 
Xc in the limit N — » oo. 

Phase space representation of projections of the ground 
state into subspaces demonstrate signals of the quan- 
tum phase transition for finite A 7 ", as shown in fig- 
ure [3] The population Husimi function Qiih^h) = 

2tt 2tt 

I I ^2^W ^l^)! 2 i s a quasi-probability distribu- 
00 

tion of the average occupations in the three local 
modes, while the function <£(</>i, </>2) = |(0i, <j)2\^)\ 2 ^ with 
|0i,02) = E e- in ^- in ^\n u n 2: n s } : repre- 

n 1 +n 2 +n 3 =N 

sents the probability distribution of collective phase dif- 
ferences [23] between the local condensates. Considering 
X = [i = in Fig. [3ja, the exact ground state is the coher- 
ent state \N;wi = W2 = 1), centered at the equilibrium 
point 1 + , showing the non-preferential occupation of the 
local modes. For increasing value of x m Fig. [3jb, we 
observe the state expansion in the population subspace, 
accompanied by squeezing in the conjugate subspace, re- 
lated thus to phase squeezing. Note that the coherent 
states are the most localized states in the phase space, 
but not necessarily in a subspace. When x — Xc there 
is no abrupt change in the representations of the ground 
state, since such changes should happen in a continuous 
way for finite N. When the self-collision parameter takes 
values slightly larger than Xc we notice a profound change 
in the ground state, shown in Fig. [3]c as a trifurcation of 
the occupational distribution, each of the components is 
more squeezed than the original coherent state. This is 
a accompanied by the emergence of an interference pat- 
tern in the conjugate subspace. This fragmentation of 
$ represents the superposition in states of preferential 
occupancy in each local mode, also visible in the behav- 
ior of Qi. Notice that the ground state representations 
in phase space display a smooth transition from RO to 
MST regime, in contrast to the abrupt transition shown 
in the macroscopic limit. 

Now we return to the discussion of particle entangle- 
ment, as given by the generalized purity, and the signal- 
ing of the QPT [3 j. In figure jij we show d ^ (3) calculated 
at the ground state as a function of x, considering sev- 
eral values of N. In the absence of bosonic collisions 
the ground state is a coherent state and consequently its 
purity value is equal to one. For x < Xc? the region re- 
lated to the expansion-compression of the ground state 
representations, we see that V su (s) decreases slowly with 
increasing \. However, for values of x slightly above 
Xc we observe a rapid decay of V su ^ caused by the 
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ground state fragmentation associated with the tunnel- 
ing suppression. The loss of purity is faster for larger 
number of trapped bosons, suggesting a scaling property 

dV su (3) 



between ^} 3) and N. We define the scalable quan- 
tum critical parameter xl(N) as the self-collision pa- 

dV su(3) 



rameter value that minimizes ^ {S) for a specific num- 
ber of particles. We observe that the minimum value 
of dV ^} 3) becomes more pronounced for increasing TV, 
while xj?(iV) moves to the left, toward Xc- The values 
of Xc(N) obtained from figure [4] suggest a linear rela- 
tionship between ln[x-?(7V) — Xc] and ln(iV). A linear in- 
terpolation of the data provides the following power-law: 



X q c (N) - xc 



O 1.2±0.2 



N' 



-0.99±0.05 



Therefore, the con- 



vergence law of Xc(N) to Xc = 2 as TV — » oo is character- 
ized. Also P sn (3) signals the QPT correctly, considering 
the scaling behavior for finite N [26] . 
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Figure 4: (color online) We show — ^ 
state as a function of x, for Q = — 1, fi 
N. Xc is indicated by the dashed vertical line 



— calculated in the ground 
= and several values of 



Although the critical bifurcation parameter x+ is dis- 
tinct from the critical transition parameter x c , the emer- 
gence of new fixed points in the semiclassical dynam- 
ics shows a clear relation with the occurrence of the 
QPT. This same characteristic has been observed in 
the two local modes model [24], but with some impor- 
tant differences. The two- mode model presents a sec- 
ond order QPT, accompanied by a pitchfork bifurcation 
with critical parameter equal to the critical parameter 
^tm _ ^tm _ i ^0 , in the absence of cross-collisions. 

However, in both models the equilibrium point that 
minimizes the semiclassical energy H can be described 
by a single angle on S 2 , considered as a function of 
X- For example, the restriction wi = W2 G R determines 
the unidimensional manifold of the phase space that con- 
tains the minimum energy fixed point. Therefore, we 
generically consider H restricted to only one dimension 
in the analysis of the QPT: U(6;x) = ^o(0) + X^i(0); 
where Ho (J~Li) is the portion of the semiclassical Hamil- 
tonian obtained from the linear (quadratic) terms of H 
in the generators of the dynamical group. The posi- 
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Figure 5: (color online) On the left the solid (blue) lines show the 
curves of Omin (above) and 7i(0mi n ) (below) considering the two 
local modes model. On the right we show the same curves for the 
three-mode model. Remark that differently from the two modes, 
there is an discontinuity in for the three modes. 



tion of the semiclassical ground state Omin = #mm(x)> 
which satisfies the necessary conditions ^\ e _ e = 

and 



d0 2 



> 0, and the minimum energy curve 

- min __ 

Uminix) = 7-L(6min(x),x) are shown in figure [5] for the 
two and three- mode models [27] . The function 9min(x) 
possesses discontinuous first derivative at x^ M for the 
two-mode model, while for the three-mode model we ob- 
serve a discontinuity in the curve of m i n (x) itself at the 
critical value Xc- The different types of bifurcation are 
responsible for the distinct behavior of the curve rnin (x) 
in both models, resulting also in different orders of tran- 
sition. The first derivative of H m in with respect to the 
parameter x is given by: 



dU r , 



dx 



on 

86 



dd n 



dx 



dU 

dx 



— Hl(@min) 



(7) 

Therefore, the continuity of 6 m in implies that dn ( ^ in is 
also continuous, since T-Li{0 m in) is a continuous function 
of Omin- In general, the reciprocal is also true, because 
H>i(Qmin) is discontinuous in the case of 9 m in also discon- 
tinuous, resulting in a first order QPT. The discontinuity 
in the derivative of m in can be identified only in the sec- 
ond derivative of %min'- 



d?l~Lmin dT~L\ 



dx 2 



dO 



dOr, 



dx 



(8) 



So, in general, the discontinuity of ££ n implies dis- 
continuity only for the second derivative of H m i n: result- 
ing in a second-order QPT. 

In conclusion we have described the dynamics of a 
BEC trapped in a symmetric triple well potential in a 
three mode approximation and show in details the way 
in that the ground state of the model undergoes a QPT. 
A time dependent variational principle using SU(3) co- 
herent state allows for a system of semiclassical equations 
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that enables one to find the fixed points of the model and 
to investigate how the lowest energy fixed points change 
as the collision parameters of the model are varied. Thus 
particle entanglement signals correctly the presence of 
QPT as one crosses distinct dynamical regimes. We have 
shown that the ground state of the model is related to 
the subdynamics of the twin condensate, a regime whose 
dynamics falls into the SU (2) group. We also show that 
the increase in the entanglement of the ground state is 
associated to phase compression (Rabi oscillation) and 
to number compression (self-trapping). In this way we 
relate the loss of coherence in the fundamental state to 
"spin squeezing". Remark that the squeezing is not the 
only effect that leads to the decrease of the generalized 
purity. Whenever there is a departure from a coherent 
state the purity decreases, but this behavior unambigu- 



ously describe particle entanglement in the system. The 
more entangled are the particles the smaller is the pu- 
rity associated to the pertinent algebra. We remark that 
this three symmetric well model allows for a better char- 
acterization of QPT than the case for double well, pro- 
viding not only the dynamical transitions but signaling 
the presence of squeezing and thus could arouse exper- 
imental investigation with existent trapping technology 
for this phenomenon. This amount of squeezing and en- 
tanglement, if controllable, could provide a new scenario 
for testing quantum physics fundamental questions. 
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